Intrinsic resistivity and the SO(5) theory of high-temperature superconductors 
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The topological structure of the order parameter in Zhang's SO(5) theory of superconductivity 
allows for an unusual type of dissipation mechanism via which current-carrying states can decay. The 
resistivity due to this mechanism, which involves orientation rather than amplitude order-parameter 
fluctuations, is calculated for the case of a thin superconducting wire. The approach is a suitably 
modified version of that pioneered by Langer and Ambegaokar for conventional superconductors. 
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Zhang's approach to the physics of the high-temp- 
erature superconductors |l| is rooted in the observa- 
tion that the phase diagram of these materials contains 
nearby regions of antiferromagnetism and superconduc- 
tivity. This has led Zhang to propose a low-energy ef- 
fective description of these materials that combines the 
U(l) symmetry sector associated with the superconduc- 
tivity with the SO (3) symmetry sector associated with 
the antiferromagnetism, SO (5) being the minimal sym- 
metry group that admits this combination. In the result- 
ing description, the state of the system is characterized, 
locally, by a five-dimensional "superspin',' subject to a 
symmetry-reducing term that favors one or other of the 
sectors, U(l) or S0(3). The superspin is constrained to 
have unit magnitude, which is appropriate for tempera- 
tures sufficiently low that fluctuations in the magnitude 
of the superspin are negligible. The possibility of ro- 
tations between the U(l) and SO (3) sectors is where a 
number of the unusual consequences of Zhang's approach 
lie |-|. 

It has long been appreciated |^J^ that in conven- 
tional superconductors topologically accessible fluctua- 
tions in the amplitude of the superconducting order pa- 
rameter provide an intrinsic mechanism via which su- 
percurrent can be dissipated in a thin wire. A uniform 
current-carrying state characterized by a specific uniform 
phase-gradient along the wire is only metastable, thermo- 
dynamically, so that by undergoing amplitude-reducing 
thermal fluctuations the system can decrease its current. 
The purpose of the present Paper is to investigate a 
related intrinsic dissipation mechanism appropriate for 
the superconducting state of the S0(5) model of high- 
temperature superconducting materials. In this case, the 
relevant dissipative process is a fluctuation in the orienta- 
tion of the superspin, during which the system temporar- 
ily becomes antiferromagnetic along a small segment of 
the wire, allowing the analogue of a phase-slip process to 
occur. Such fluctuations are equivalent to the passing of 
superconducting vortices with antiferromagnetic cores ||] 
across the wire. Related dissipative processes have been 
addressed in the context of superfluid "^He-A and (for 
appropriate values of the gradient coupling constants) 
thin tubes of nematic liquid crystal [B| . By constructing a 



version of the approach to intrinsic dissipation pioneered 
by Langer and Ambegaokar for conventional supercon- 
ductors 1^, suitably modified for the SO (5) model of su- 
perconductivity, we shall estimate the free-energy barrier 
for this type of fluctuation and, hence, arrive at an esti- 
mate of the current-voltage relationship for a sufficiently 
thin wire. 

In Zhang's theory of high-temperature superconduc- 
tivity and antiferromagnetism, then, the local state of 
the system at the spatial position r is determined by a 
five-component unit-vector field n(r) such that the com- 
ponents n^, and rr' together specify the antiferromag- 
netic Neel vector, and the components n'^ and it' together 
specify the amplitude and phase of the superconducting 
order. The constraint on the magnitude of n implements 
the notion that the enhancement of antiferromagnetic 
order is necessarily accompanied by the diminution of 
superconducting order, and vice versa. The free-energy 
density / comprises an isotropic gradient term along with 
a symmetry-reducing term: 
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Here, the (real-space) subscript /i runs from 1 to 3 
(or X, y, z), repeated indices being summed over. 
The (chemical-potential dependent) parameter g governs 
whether the stable homogeneous state is superconducting 
{g < 0), as we select here, or antiferromagnetic (g > 0). 
We parametrize n via the angles {9, 0, ip, x} such that 



— sin 9 cos cos x, 
~ sin 9 cos tp sin x, 
tV^ = sin sin ^/;, 
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The angle 9 measures the relative amount of antiferro- 
magnetic versus superconducting order (without regard 
for the orientation of the antiferromagnetism or the phase 
of the superconductivity), with = in the purely su- 
perconducting state. Furthermore, <j) is the phase of the 
superconductivity, and '0 and x specify the orientation 
of the antiferromagnetic Neel vector. 
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Let us suppose that the system consists of a wire of 
length L, sufficiently long and narrow that we may as- 
sume that {9, (j), ip, x} do not vary in the directions trans- 
verse to the wire (i.e. in the y or z directions). Under 
these circumstances / becomes 



f = ip/2){{djf + (d. 



)2 cos^^ 



(9,^)2 sin2 
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(2a) 



(|g|/2) {{dr9f + {dr(l>f cos^9 + {dri^f Slu'l 
+ {drxf sin26' cos^?/. + sin^e*} . 



(2b) 



To obtain Eq. (gb|) we have exchanged the independent 
variable x for its dimensionless counterpart r, defined via 
r = x/^TT, where = \/pl\g\ is the correlation length 
for antiferromagnetic fluctuations, which sets the length- 
scale for dissipative events. Furthermore, we have defined 
the function 9 such that 9(t) = 9{x), and similarly for (p, 
ij) and X- The quantity i will denote the dimensionless 
length of the wire, i.e. £ = L/^^. 

In order to estimate the rate at which current-dis- 
sipating processes occur we compute the height of the 
free-energy barrier opposing them. To do this we follow 
Langer and Ambegaokar and seek the metastable (i.e. 
uniform, current-carrying) states between which the sys- 
tem fluctuates, and the transition (i.e. unstable saddle- 
point) states through which the system passes as current 
is dissipated. Both classes of states, metastable and tran- 
sition, are stationary configurations of the free energy, 
and therefore obey the corresponding Euler-Lagrange 
equations. These equations may be simplified, how- 
ever, due to the homogeneity (i.e. r-independence) and 
gauge-invariance in the superconducting sector (i.e. 0- 
independence) of /. The former symmetry leads to the 
existence of the first integral 

i)"^ + 4^ cos^e + ^'^ sm^9+x^ sm^e cos^i:-s\r?9 = e, (3a) 

where overdots denote derivatives with respect to r, and 
e is constant. The latter symmetry leads to the existence 
of a cyclic coordinate, so that 



icos^e = I . 



(3b) 



Here, the dimensionless supercurrent density /, in terms 
of which the dimensionful supercurrent density J is given 
by J = 2ep//?i^7r, is constant. By replacing in Eq. ( |3a| ) 
with /, using Eq. (3b), we obtain 



9^ + 



cos'^9 



+ Tp^ s\v?9 + x^ sin^e cos^?/) - sin^e* = e . (4) 



The relevant supercurrent-carrying metastable states 
are solutions of Eqs. (3) for which 9{t) = (and there- 
fore Tp and X pl&y no role) and (p — I. We refer to 
these states as "uniformly winding" because in them <p 
increases linearly with position along the sample. As 
for the transition states, we assume that they possess 
appreciable antiferromagnetic order only within a small 



segment of the wire, and that elsewhere 9 is negligibly 
small. Furthermore, the orientation of the antiferromag- 
netic order generated during the fluctuation is uniform, 
because transition states with inhomogeneous antiferro- 
magnetic orientation would have higher free energy, and 
would correspondingly occur less frequently. Therefore 
we only consider transition states for which i/) = x = 0. 
Thus, by applying Eq. (^ far from the antiferromagnetic 
region we determine that 



/2 = e. 



(5) 



Let us denote the extreme value of 9{t) by 6*0, and let 
us suppose that it occurs at the location r = tq. At this 
location ^ = 0, so that Eq. (^), combined with Eq. (||), 
gives 



9^ ^ sm^9 - iim^L 
I? = cos26'o . 



By integrating Eq. (pa) we find 



sin 6'(r) = sin 9o y I ~ tanh [(r — tq) sin 9q] 



(6a) 
(6b) 



(7) 



Finally, by using Eq. (3b), along with Eq. (|7|), we obtain 
0(r) - 0(to) = /t / dr'/ cos'' 9{t') (8) 



= /t (r — To ) + arctan {tan tanh [(r — tq) sin 9q] } , 




FIG. 1. Parametric plot of a typical (partially antiferro- 
magnetic) transition state (full curve) on the surface of the 
unit 2-sphere in the space spanned by {n*, n^, sin 6}. Trajec- 
tories of the (purely superconducting) metastable states lie 
on the equator of this sphere, for which ^ = 0. 

Thus we see that far from tq (i.e. far from the center 
of the fluctuation) the transition states wind uniformly, 
as do the metastable states. However, in the transition 
states the order parameter undergoes an orientational 
distortion over a length of order that takes 9 out of the 
9 = plane. Thus, a given transition state possesses a 
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region of antiferromagnetic order in a fixed but arbitrary 
direction. Indeed, there is a family of symmetry-related 
transition states generated by changing this direction, 
just as there is a family of transition states generated by 
the (arbitrary) location of tq and the (arbitrary) over- 
all phase. An example of a transition state is shown 
in Fig. |l|. Hence we see that the relevant fluctuation 
for SO (5) superconductors is one in which supercurrent 
can be dissipated via a thermally activated process in 
which a "loop" of order parameter passes over the order- 
parameter sphere, so that the total phase difference (and 
hence supercurrent in the resulting metastable state) is 
reduced by 27r. Intuitively, it seems reasonable that if 
the "anisotropy" g is not too large (it being adjustable 
by changing the chemical potential) then this type of 
fluctuation should be less energetically costly than ampli- 
tude fluctuations, and should hence provide the dominant 
pathway for dissipation. Having determined the form 
of the metastable and transition states, we now identify 
the particular transition state through which the system 
passes as current is dissipated from a given metastable 
state. In passing from a metastable state to the cor- 
responding transition state, en route to diminishing the 
current, the superspin nowhere becomes entirely antifer- 
romagnetic (i.e. \9\ < 7r/2). [The transition state rep- 
resents a configuration in which a loop is about to pass 
over the order parameter sphere, but has not yet done 
so.] Thus, if we consider a transition from one metastable 
state (which we refer to as the "upper metastable state" ) 
to a metastable state with 27r less total phase difference 
(the "lower metastable state" ) then the actual loss of the 
phase will occur in passing from the transition state to 
the lower metastable state. Therefore, the total phase- 
difference across the sample A(j) is the same in the tran- 
sition state and the upper metastable state, so that 



dr ( 



A0: 



dr ( 



(9) 



where (jj^^ (t) is the phase of the uniform metastable state 
and (f>t{T) is the phase of the transition state. This for- 
mula provides a connection between the current in the 
upper metastable state and that in the transition state 
It- The left- most term in this equation is given by I^l, 
which follows readily from the uniformly twisting char- 
acter of the metastable states. The right-most term may 
be evaluated via Eg. (||), which gives It£ + 26o, or equiv- 
alently, using Eq. (6b), It£ + 2arccos/t. Thus we arrive 
at an implicit equation for It in terms of /+: 



I+e-It£^ 2arccos/t > 0. 



(10) 



Next we use Arrhenius rate-law considerations to de- 
velop the transition rate. To do this, we need expres- 
sions for the free energies of the upper (-I-) and lower (— ) 
metastable states, as well as of the transition state con- 
necting them. From Eq. ( p^ ) and the form of the upper 
and lower metastable states we integrate over the volume 



of the wire to obtain expressions for the free energies 
of the these states: 



F± = (|5|/2)yl^.£(/±)'. 



(11) 



Here, (= /+ — 2Tr/£) is the current in the lower 
metastable state and A is the cross-sectional area of the 
wire. Similarly, the free-energy density of the transition 
states ft may be obtained using Eqs. (pb|), ( ^ ) and (|6^): 

/t = (|ff|/2)(/2 + 2sin20). (12) 

By integrating Eq. ( p^ over the volume of the wire we 
arrive at an expression for the free energy Ft of the tran- 
sition state: 

Ft^{\g\/2)A^^{£I^ + 2 dr sin^ 6 (t)} (13a) 
= i\g\/2)A^^[£I^ + A,/r^}. (13b) 



By using Eqs. (|llj) and (13b) we obtain an expression for 
the free energy barrier AF for current dissipation. 



AF = Ft-F+ = 



\9 



■A^^{£I^-£{I+r+4.Jl^?}, 



which may be simplified by using the relation (jl^) be- 
tween It and I^. By restricting our attention to states 
in which the phase winds many times along the wire (i.e. 
2 arccos/t <C It £) we obtain I"^ — (/^)^ — — 4/t arccos/t, 
and thus AF becomes 



AF 



2 It arccos L 



+ 2^1^}. (14) 



It should be noted that /t = 1 is the critical current, in 
the sense that uniformly twisted states with larger values 
of the current are unstable rather than metastable. For 
later use, we note that for currents slightly smaller than 
critical AF in the SO (5) model is given by 



AFc^2^\g\A^^{l~Itf'\ 



(15) 



whereas the corresponding expression for conventional 
superconductors [p| is 



AF ^ V2 (8/3)5/4 AC _ _ ^)5/4 ^ 



(16) 



Here, ((?„ — ffs) represents the free-energy cost of ampli- 
tude fluctuations, and is the analogue of the parameter g 
within the present theory, and is the superconducting 
fluctuation correlation length. 

To compute the rates of current-decreasing and 
current-increasing fluctuations we follow LA by assuming 
that these rates depend exponentially on the free-energy 
barrier heights. Specifically, the rate r(/+ I~) at 
which current-decreasing fluctuations occur is given by 



r(/+^/,;) = roexp{-/3AF) 



(17a) 



Fo exp <^ -p \g\ A^^ ~2It arccos/t + 2Jl - I^ 
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where (3 = l/k^T measures the inverse temperature 
and To is an attempt frequency for dissipative fluctua- 
tions. Current-increasing fluctuations, in which the sys- 
tem passes from the lower metastable state to the up- 
per metastable state, also occur. However, these fluctu- 
ations have a higher barrier opposing them (i.e. for them 
AF AF + F^ — F~), and thus occur at a lower rate. 
By using Eq. (^l|) we therefore obtain 



(17b) 



Toexjp<-(3\g\A^^[-2It arccos It+2Jl-I^+ 2ttI, 



The application of a voltage difference AV between the 
ends of the sample causes to increase linearly with 
time [i.e. d{A(j>)/dt = Aire AV/h, where h is Planck's 
constant and — e is the electronic charge] In steady 
state, this will be balanced by the net rate Fnot at which 
current-decreasing fluctuations occur, i.e., 



= r(/+ 



-r(/^ 



/-) = 



1 d{A(f)) 
2tt dt 



(18) 



The difference between It and is of order , and 
we shall henceforth neglect it; thus we simply refer to 
the current /. This leads to the following expression for 
the voltage between the ends of the wire in terms of the 
current along it: 



Al/= (Ve)ro smh{nl3\g\A^^I} 

X exp|-/3|5|^^^ (^-2/arccos/ + 2v/l - P 



(19) 
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The corresponding expression for conventional supercon- 
ductors may be found in Eq. (2.12) of Rcf. As the 
length-scale of the fluctuation is it is natural for AV to 
depend exponentially on /3|g|^^7r (i.e. the condensation 
energy in a correlation length per kBT). Following Mc- 
Cumber and Halperin |l^] , we estimate Tq to be of order 
N/t, where N = AttL/^^-^ reflects the range of possible 
locations along the wire at which the fluctuation may oc- 
cur (i/^^), as well as the variety of possible Neel- vector 
orientations (47r). For conventional superconductors r 
diverges as T ^ Tc- For the SO (5) model we expect r to 
diverge as the chemical potential approaches its critical 
value. 

It may be useful to compare expression (^9|) for Ay 
with that obtained in Refs. ||^,|l^ for the case of con- 
ventional superconductors. Within the SO (5) model the 
expression for AV in terms of the dimensionful supercur- 
rent density J is given, for J ~ J^. , by 

AV ~ (h/e) Tq sinh{/3A/iJ/4e} (20) 
X exp {-f3A ({hJ/Ae) + 2^ \g\ (1 - J/ Jc)i ) } , 

whereas the corresponding expression obtained in 
Refs. [pp^ is given, for J ~ Jc, by 



AV ~ {h/e) VLq smh{PAhJ/Ae} (21) 
5{-/3A((/iJ/4e) + V2(8/3)4c(<7n-.9s)(l-J/Jc)*)}. 



xexp-^ 



(We have ignored the current dependence of the attempt 
frequency f^o obtained in Ref. ||l^].) We see that in this 
regime the primary qualitative distinction between for- 
mulas (|20|) and (gl|) comes from the distinction between 
the nature of the fluctuations, and the resulting differ- 
ence between the barrier heights. 

We conclude by noting that for conventional supercon- 
ductors LA 1^ have obtained, as a condition for the va- 
lidity of their approach, the constraint that the wire be 
thinner than the temperature-dependent superconduct- 
ing correlation length ^. (To accomplish this they con- 
sider the second variation of the free energy, and require 
that the transition state be unstable in only one direction 
in configuration space.) We emphasize that the corre- 
sponding condition in the context of the S0(5) theory of 
high-temperature superconductivity is that the wire be 
thinner than the correlation length for orientation fluc- 
tuations of the superspin, i.e., ^jr- This length is ex- 
pected to diverge when the chemical potential approaches 
the superconductor-to-antiferromagnet phase boundary. 
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